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ABSTRACT 


This  report  describes  a  program  of  work  on  beam/plasma  interaction. 
Both  electrostatic  and  electromagnetic  wave  amplifying  mechanisms  are 
under  investigation.  For  the  former,  studies  in  the  absence  of  a  static 
magnetic  field  are  directed  towards  verifying  the  theory  for  the  cases 
of  finite  beam/infinite  plasma,  and  beam/surface  wave  amplification, when 
transverse  modulation  is  applied,  A  dipole  resonant  coupling  system  for 
such  interactions  is  under  study.  Two  distinctly  different  lines  are 
being  followed  for  interactions  in  tl.  'esence  of  a  static  magnetic 
field;  Electrostatic  cyclotron  harm  tic  wave  interaction  is  being  ex¬ 
amined,  both  theoretically  and  experimentally,  and  the  potentialities  of 
electromagnetic  wave  growth  in  the  "whistler  mode  are  being  investigated. 


in 


FOREWORD 


This  contract  represents  a  three-year  program  of  research  on  "Fast 
Wave  Beam/plasma  Interactions"  which  is  proceeding  in  the  Institute  for 
Plasma  Research,  Stanford  University,  under  the  direction  of  Prof.  F.  W. 
Crawford  as  Principal  Investigator.  The  work  is  part  of  PROJECT  DEFENDER 
and  was  made  possible  by  the  support  of  the  Advanced  Research  Projects 
Agency  under  Order  No.  695*  If  Is  conducted  under  the  technical  guid¬ 
ance  of  the  U.  S.  Army  Electronics  Command.  This  is  the  fifth  Quarterly 
Report,  and  covers  the  period  1  March  to  31  May,  1967. 
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I .  INTROTUCTION 


The  wave  amplification  effect  associated  with  the  interaction  of 
an  electron  beam  and  a  plasma  has  attracted  considerable  attention  over 
the  last  few  years,  particularly  from  microwave  tube  specialists  to 
whom  such  interactions  offer  possibilities  of  constructing  very  high 
gain  devices  which  should  be  electronically  tunable  over  wide  frequency 
ranges.  Since  the  plasma  plays  the  role  of  a  conventional  slow-wave 
structure,  the  interaction  region  should  be  free  of  metallic  structures, 
a  particularly  significant  characteristic  if  millimeter  wave  operation 
is  envisaged. 

The  work  being  carried  out  under  this  contract  is  directed  towards 
utilizing  the  beam/plasma  amplification  mechanism  in  microwave  device 
applications.  So  far,  despite  the  efforts  of  many  groups,  it  has  not 
been  found  possiblo  to  realize  this  potential  fully.  The  most  serious 
obstacles  to  progross  are  that  efficient  coupling  of  an  rf  signal  into 
and  out  of  the  interaction  rogion  has  been  found  difficult  to  achieve, 
and  that  tho  amplifiers  arc  frequently  very  noisy  compared  with  most 
conventional  microwave  tubes.  The  necessity  of  providing  the  means  of 
plasmn  g  ntion  within  tho  device,  and  the  presence  of  a  relatively 
high  background  gas  prossuro,  add  constructional  problems  beyond  those 
normally  encountered  with  vacuum  tubes.  Although  satisfactory  engineer¬ 
ing  solutions  to  thoso  latter  difficulties  could  certainly  be  found, 
the  coupling  and  noiso  problems  still  require  considerable  further  study 
to  dotermino  whethor  competitive  devices  can  be  developed. 

Of  the  many  Widely  differing  aspects  of  beam/plasma  interaction, 
three  have  been  chosen  for  closo  examination  under  this  contract.  The 
first  of  these  is  the  interaction  of  an  olectron  beam  with  a  plasma 
when  the  modulating  fields,  and  tho  wave  growth,  are  in  the  first  azi- 
muthally-varying  mode.  Since  with  transverse  modulation  several  inter¬ 
esting  interaction  and  coupling  mechanisms  become  possible,  it  is 
intended  that  a  major  portion  of  the  work  under  this  contract  should 
consist  of  a  thorough  investigation  of  such  phenomena. 

Most  previous  work  has  been  concerned  with  the  theoretical  descrip¬ 
tion  and  demonstration  of  beam/plasma  interaction  mechanisms  that  can 


be  derived  from  cold  plasma  theory,  i.e,,  from  theory  in  which  it  is 
assumed  that  the  plasma  electrons  have  no  thermal  or  directed  motions, 
and  that  the  irjected  beam  is  monoenergetic.  When  a  oc  magnetic  field 
is  present,  microscopic  theory,  in  which  single-particle  behavior  is 
followed,  predicts  a  much  wider  range  of  amplification  mechanisms.  Some 
of  these  are  simply  modifications  of  those  occurring  in  the  absence  of 
the  magnetic  field,  while  others  involve  interaction  of  beams  with  trans¬ 
verse  energy  with  slow  "cyclotron  harmonic  waves."  This  constitutes  our 
second  area  of  interest,  i.e.,  that  of  wave  growth  in  magnetoplasmas 
3n  the  electron  beam  has  a  substantial  component  of  transverse  energy. 

Our  third  area  of  interest  is  in  electromagnetic  wave  amplification. 
Theoretical  studies  show  that,  in  addition  to  electrostatic  wave  growth 
phenomena  such  as  those  just  described,  there  is  the  possibility  of 
obtaining  appreciable  growth  in  the  "whistler"  mode  when  an  electron 
beam  with  transverse  energy  interacts  with  a  magnetoplasma.  This  mode 
is  a  right-hand,  circularly-polarized  electromagnetic  wave,  i.e.,  its 
electric  field  vector  rotates  in  the  right-hand  sense,  which  is  also 
( conventionally)  the  sense  of  rotation  of  the  electrons  about  the  mag¬ 
netic  field  lines.  If  a  beam  with  transverse  energy  is  moving  aiong 
the  field  lines,  there  is  consequently  a  possibility  of  energy  being 
transferred  from  the  electrons  to  the  wave,  and  hence,  for  wave  ampli¬ 
fication  to  occur.  The  purpose  of  our  work  is  to  demonstrate  this  type 
of  interaction,  and  to  examine  its  potentiality  for  coupling  to  slow- 
and  fast-wave  circuits.  Here  fast-wave  is  interpreted  to  mean  that 
the  phase  velocity  of  the  wave  is  of  the  order  of  the  velocity  of  light. 

Previous  quarterly  reports  (QR)  have  described  the  backgrr  ad  for 
each  of  the  topics  in  detail.  Progress  made  during  the  reporting  period 
will  be  described  in  the  succeeding  sections. 
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II.  BEAM/PLASMA  AMPLIFICATION  WITH  TRANSVERSE  MODULATION 

In  previous  QR, .  various  general  types  of  electrostatic  beam/ 

Pi  asma  interactions  have  been  discussed  in  detail.  The  particular  con¬ 
ditions  of  interest  to  us  under  this  project  are  those  of  a  finite  beam 
interacting  with  a  finite  plasma  surrounded  by  one  or  more  dielectrics 
and  a  metallic  waveguide.  We  may  distinguish  two  cases:  one  in  which 
the  beam  diameter  is  small  compared  with  the  plasma  diameter,  and  one 
in  which  the  beam  fills  the  plasma.  In  the  latter  case,  the  interaction 
occurs  between  the  beam  and  surface  waves  that  may  propagate  along  the 
column.  Our  aim  is  to  study  both  types  of  interaction,  both  theoretic¬ 
ally  and  experimentally,  for  modulation  applied  transverse  to  the  beam. 
This  mode  of  interaction  has  inherently  interesting  features  in  relation 
to  improved  coupling  between  the  interaction  region,  and  structures  out¬ 
side. 

In  QR  1+,1  the  question  of  the  optimum  choice  of  plasma  to  waveguide 
radius  ratio,  and  the  maximum  beam  voltage  for  controlled  mode  selac’ion 
(m  =  0  or  l)  were  discussed,  together  with  preliminary  experimental  re¬ 
sults  on  input/output  coupling,  propagation  of  the  amplified  signal,  and 
the  appearance  of  ion  oscillations.  These  studies  have  been  pursued 
during  the  current  reporting  period. 

(A)  Theoretical  Studies  of  Beam/plasma  Interaction. 

For  the  two  beam/plasma  radius  ratio  limits  to  be  studied,  it 
is  expected  that  when  the  ratio  is  small,  the  spatial  growth  rate  will 
be  high,  but  that  the  electric  field  near  the  glass  wall  of  the  tube  and 
in  the  waveguide  will  be  weak,  making  strong  coupling  difficult  to  real¬ 
ize.  The  other  limit,  with  the  beam/plasma  radius  approaching  unity, 
provides  much  stronger  electric  field  at  the  dielectric  interface,  and 
therefore  better  coupling.  Coupling  Q's  up  to  50  have  been  observed  in 
this  configuration  and  reported  in  previous  QR.  Due  to  the  larger  inter¬ 
action  region,  i.e.  lower  beam/plasma  frequency  for  the  same  current,  the 
gain  is  somewuat  lower.  Amplification  in  this  mode  may  also  be  more 
strongly  affected  by  radial  plasma  inhomogeneity. 
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In  the  quasistatic  approximation,  the  relevant  dispersion  relation 
for  the  geometry  of  Fig.  1  is 
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In  order  to  mnko  a  detailed  investigat ion  of  Eq ,  (5),  a  program  was 
first  written  to  find  the  renl  roots  of  the  dispersion  relation.  These 
are  shown  in  Fig.  1.  It  will  be  seen  that  for  the  in  ■  0  mode,  the 
fast  surface  wave  in  the  beatn/plasma  interaction  case  approaches  the 
plasma  dispersion  for  the  long  wavelengths,  whereas  for  short  wavelengths 
it  has  for  asymptote  the  fast  beam  surface  wave.  For  (“^q)  <  1  there 
are  two  pure  real  roots,  the  second  having  for  asymptotes  the  slow  beam 
surface  wave  and  the  plasma  surface  wave. 

Our  current  efforts  are  directed  to  determiring  the  complex  roots 
of  the  dispersion  relation  for  both  the  m  ■  0  and  m  ■  1  modes. 

(B)  Experimental.- Studies  of  Beam/plosma  interaction. 

In  QR  4,  computations  of  tube  geometry  such  that  the  beam 
characteristic  intersects  the  m  =  1  plasma  mode  only  were  presented, 
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and  It  waa  pointed  out  that  thin  could  be  assured  in  two  way*:  either 
by  increasing  the  bean  voltage,  or  by  reducing  the  ratio  (c/»)  (nee 
Fig.  1).  During  the  quarter,  a  new  boom/waveguide  system  has  been  con¬ 
structed  with  (c/e)  •  l.y>  ,  end  designed  for  interaction  with  the 
n  -  1  code  only  for  Vb  a  <6$  V.  This  is  shown  in  Fig.  2.  Facilities 
have  been  provided  for  excitation  of  the  a  ■  0  and  a  ■  1  modes  sep¬ 
arately,  and  the  input  coupler  can  bo  rotated  to  distinguish  the  node 
of  excitation  actually  occurring  under  given  conditions.  One  end  of 
the  tube  has  been  aade  non-reflecting  by  provision  of  absorbing  mater¬ 
ial  and  a  coating  of  Aquadag'.  The  first  swasuremcnls  with  the  new 
systca  concerned  the  dc  stability  of  the  discharge  and  choico  of  the 
best  operating  point.  It  was  found  that  evon  with  a  two-grid  system, for 
control  and  acceleration  of  the  beam,  the  discharge  could  bo  space- 
charge  Halted  for  only  n  short  tiac,  and  that  long  tera  stability  of 
the  column  could  bo  achieved  through  teaperature-liaited  operation  only. 
The  best  DC  operating  point  was  found  to  bo  Vb  ■  1^0  V  ,  and 
0b  *  *p)  ■  1«*  ,  tor  which  the  beaa  line  intersects  both  the  m  ■  0 
and  a  -  1  modes,  however.  To  ensure  stability  of  the  plasma  density, 
a  highly-stoblc  DC  heater  supply  is  used. 

Since  the  best  operating  point  still  allowed  propagation  and  growth 
in  the  a  ■  0  mode,  a  number  of  neasureaents  were  made  of  its  charac¬ 
teristics.  It  is  essential  to  know  these,  and  to  avoid  such  effects  os 
beaa  break-up  due  to  noise  in  the  ra  ■  0  mode  while  studying  growth  of 
modulated  signals  in  the  ra  ■  1  mode.  Some  representative  data  on  amp¬ 
litude  and  phase  Blong  the  tube  i.-c  shown  in  Fig.  5.  It  can  be  seen 
that  the  growth  rate  is  quite  low;  less  than  1  dB/cm.  When  the  compu¬ 
tations  of  the  coupled  roots  of  the  dispersion  relation  are  completed, 
wo  will  be  able  to  check  the  agreement  between  simplified  theory  and 
those  experimental  results. 

The  beam  break-up  phenomenon  referred  to  manifests  itself  as  shown 
in  Fig.  k,  by  a  sudden  change  of  the  amplitude  of  the  spatiolly  growing 
signal,  ind  somewhat  increased  phase  shift. 

When  the  break-up  is  caused  by  the  modulating  signal,  the  bicak-up 
region  moves  towards  the  cathode  as  the  plasma  frequency  is  approached 
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III.  ELECTROSTATIC  WAVE  AMPLIFICATION  IN  MAGNETOPLASMAS 


Whon  the  beam  and/or  plasma  have  directed  or  thermal  motions  in  the 
tronsvorse  and  axial  directions,  it  is  necessary  to  derive  the  appropri¬ 
ate  dispersion  relations  using  a  Boltzmann  equation  formalism.  The  re¬ 
sults  of  doing  so  were  discussed  rather  generally  in  QR  1  where  it  was 
polntod  out  that,  for  a  high  enough  value  of  the  parameter  (o>  /u>  )  i.e. 
the  ratio  of  beam  plasma  frequency  to  electron  cyclotron  frequency,  even 
an  ion-neutralized  electron  beam  could  be  unstable,  and  that  in  the  pres¬ 
ence  of  a  backpround  plasma  the  instability  threshold  for  the  beam  density 
could  bo  reduced.  The  purpose  of  this  project  is  to  investigate  such 
interactions,  and  to  determine  their  potentialities  for  microwave  appli¬ 
cations. 

Numerous  theoretical  predictions  of  the  instabilities  have  been 
mdo  at  Stanford  and  olsewhere.  Basically,  the  theory  predicts  growth 
in  passbands  centered  on  the  electron  cy  lotron  harmonic  frequencies 
(K»c)  *  No  furthor  computations  will  be  carried  out  under  this  project 
until  our  experimental  parameters  have  been  measured.  Those  computa¬ 
tions  carried  out  to  date  are  being  summarized  in  a  Ph.D.  thesis  being 
written  by  J.  A.  Tataronis. 

So  far,  few  controlled  experiments  have  been  carried  out  to  check 
the  theory,  though  observations  of  strong  noise  emissions  from  magneto- 
plasmas  containing  charged  particles  with  appreciable  transverse  veloci¬ 
ties  provide  significant  support  for  the  existence  of  the  predicted 
attchnnisas.  The  studies  planned  under  this  contract  are  intended  to 
provide  results  undor  refined  experimental  conditions,  and  to  put  the 
theory  on  a  firm  quantitative  basis.  In  particular,  we  wish  to  verify 
the  dispersion  relation  for  the  realistic  case  of  a  delta-function  beam 
interacting  with  a  warm  plasma. 

(A)  Experimental  Studies 

The  aim  of  the  experimental  work  under  this  project  is  to  ex¬ 
cite  growing  waves  by  means  of  an  electron  beam  injected  into  the  plasma, 
and  to  study  the  variation  of  the  growth  rate  as  i  function  of  the  lon¬ 
gitudinal  and  tronsverso  onergies  of  the  beam.  The  first,  and  simplest, 
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waj  of  imparting  transverse  energy  to  the  beam  is  to  inject  it  through 
an  increasing  magnetic  field  into  the  interaction  region.  This  does  not 
create  the  delta-function  transverse  velocity  distribution  which  would 
be  most  desirable  for  checking  against  theory.  A  more  satisfactory 
approach  is  the  use  of  a  "corkscrew"  injection  system.'  A  third  method 
which  it  was  hoped  to  apply  because  of  its  greater  flexibility  is  to 
Impart  transverse  energy  to  the  electrons  by  cyclotron  heating  in  a 
small  rf  cavity  through  which  the  beam  passes  before  entering  the  plasma 
region.  This  has  been  found  extremely  difficult  to  realize  fci  our  ex¬ 
perimental  conditions,  however.  For  convenience  in  our  initial  studies, 
the  first  method  has  finally  been  adopted. 

The  experimental  set-up  is  as  shown  in  Fig.  5.  Eight  large  coils 
have  been  spaced  so  as  to  give  an  axial  magnetic  field  uniform  to  better 
than  0.5^  over  a  legion  of  five  to  six  centimeters  centered  at  the  probe  . 
To  produce  a  rapidly  varying  magnetic  field  close  to  the  cathode,  two 
auxiliary  coils  have  been  provided.  These  reduce  the  field  in  the  vicin¬ 
ity  of  the  cathode,  but  do  not  disturb  the  uniformity  of  the  field  in 

the  region  near  the  probe  by  more  than  two  or  three  percent.  The  fields 

due  to  the  two  coil  systems  are  controlled  by  separate  power  supplies  so 

that  inhomogeneity  can  be  varied  relative  to  the  background  homogeneous 
field. 

The  theoretical  predictions  for  this  type  of  instability  indicate 
that  it  will  generally  be  absolute,  i.e.  signals  will  grow  from  noise, 
and  external  modulation  will  not  be  required  to  obtain  an  output.  This 
has  been  confirmed  qualitatively  by  the  previous  experimental  work  which 
it  is  hoped  to  refine  quantitatively  in  these  studies.  For  example  in 
a  PIG  discharge  Landauer  observed  radiation  out  to  about  the  45th  har¬ 
monic,  fitting  the  relation  (b)/0)  )  =  n  to  better  than  24.  Bekefi  and 

4  c 

Hooper  observed  strong  cyclotron  harmonic  radiation  from  a  beam  gener¬ 
ated  discharge  in  mercury-vapor.  As  here,  they  produced  the  necessary 
transverse  energy  by  magnetic  field  inhomogeneity.  Ikegimi  and  Crawford5 
have  also  made  measurements  of  radiation  near  cyclotron  harmonics  for  a 
beam-generated  mercury-vapor  plasma  in  a  magnetic  field  produced  by 
Helmholtz  coils.  The  mechanism  for  producing  transverse  beam  energy 
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FIG.  5.  Experimental  set-up  for  study  of  cyclotron  harmonic 
wave  amplification. 


was  again  an  inhomogeneous  magnetic  field 

During  the  reporting  period,  we  have  been  studying  radiation  from 
an  argon  positive  column  discharge  immersed  in  the  inhomogeneous  mag¬ 
netic  field  produced  by  the  set-up  of  Fig.  5*  Use  of  the  broad-band 
spectrum  analyzer  (HP  851)  greatly  simplifies  the  radiation  measurements 
since  it  eliminates  the  necessity  of  sweeping  the  magnetic  field. 

Typical  data  are  shown  in  Fig,  6(a)  for  the  signal  received  at  the  probe. 
The  cyclotron  frequency  is  530  MHz  at  the  probe  and  ~  5  MHz  at  the  cath¬ 
ode  surface.  Certain  peaks  in  the  spectrum  varied  in  frequency  with 
changes  in  the  discharge  current.  For  example,  analysis  of  the  varia¬ 
tion  of  Peak  A  shows  that  it  varies  with  the  square  root  of  the  dis¬ 
charge  current  and  corresponds  to  the  electron  plasma  frequency  (see 
Fig.  6(b)  ).  It  is  probably  excited  by  beam/plasma  interaction  in  an 
axisymmetric  mode.  Peaks  B  and  C  follow  more  complicated  variations 
which  are  currently  being  elucidated. 
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FIG.  6.  Noise  radiation  from  a  magnetoplasma: 

(a)  Spectral  distribution  of  the  noise. 


950 


-17- 


FIG.  6.  Noise 


IV.  ELECTROMAGNETIC  WAVE  AMPLIFICATION  IN  MAGNETOPLASMAS 


In  the  absence  of  a  static  magnetic  field,  interaction  of  an  elec¬ 
tron  beam  with  a  plasma  leads  only  to  electrostatic  beam/plasma  inter¬ 
actions  of  the  types  described  in  Section  II.  When  a  static  magnetic 
field  is  present,  there  are  additional  possibilities  of  electromagnetic 
wave  interaction.  That  of  special  importance  under  the  present  contract 
is  the  interaction  with  the  right-hand  polarized  electromagnetic  wave 
known  in  ionosphere  terminology  as  the  "whistler"  mode.  It  has  been 
demonstrated  theoretically  that  under  conditions  where  a  beam  with  trans¬ 
verse  energy  interacts  with  the  plasma,  wave  growth  in  this  mode  should 
be  possible,  and  that  experimental  situations  in  which  this  dominates 
over  the  electrostatic  growth  mechanisms  occurring  at  the  same  time 
appear  to  be  realizable.^ 

Comparatively  little  experimental  work  has  been  reported  so  far  on 
propagation  of  the  whistler  mode  in  laboratory  plasmas,  and  none  of 
this  seems  to  have  been  directed  towards  observation  of  wave  growth  due 
to  interaction  with  a  gyrating  electron  stream.  Such  a  demonstration 
forms  the  primary  object  of  this  project.  If  growth  in  the  whistler 
mode  could  be  demonstrated,  and  utilized,  it  would  offer  very  attractive 
practical  features.  In  particular,  coupling  should  be  facilitated, 
since  the  amplification  occurs  in  an  electromagnetic  mode,  i.e.,  without 
conversion  to  an  electrostatic  mode. 

The  aims  of  the  present  project  are  as  follows:  First,  to  elucidate 
the  theory  of  the  whistler-type  instabilities  in  the  simplest  geometry, 
and  then  to  extend  this  to  more  realistic  physical  conditions,  and 
second  to  demonstrate  directly  by  experiment  that  growth  can  occur  in 
this  mode. 

(A)  Theoretical  Studies. 

In  QR  4,  a  set  of  numerical  solutions  for  the  simultaneous 
equations  D(tu,k  )  =  0  ,  and  dD(o),k  )/dk  )  =  0  were  presented  with 

"  II  *' 

the  beam  velocity,  v  ,  and  the  transverse  speed  of  the  beam  particles, 
v0i  as  Parameters.  The  limiting  values  as  v^  approaches  zero  have 
been  investigated  analytically  during  this  quarter,  with  the  conclusion 
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the  brunch-point  locus  should  appro. eh  continuously  the  values  cor¬ 


responding  to  v  =  0  ,  that  is 


03  V 

03  =  03  ±  i  — ,  Pi. 


S72 

2  c 


(6) 


»»  v0|  decreases.  This  continuity  Is  shown  In  Fig.  7.  The  case  ln 

.his  limit  was  studied  by  Sudan.'  „,so  concluded  that  an  absolute 
instability  occurs.  For  a  practical  laboratory  bcam/plasma  Interaction 
the  beam  density  .11,  be  -all.  The  locus  ol  branch  points  In  the  cos,-’ 
plex-c.  plane,  as  decreases,  Is  shown  ln  Fig.  8.  The  corres¬ 

ponding  saddle-point  plot  1,  given  in  Flg.  9.  For  most  of  the  range 
studied,  a,  the  beam  density  decreases,  the  t  -mpor.r  growth  rate  In- 
creasos,  while  the  spatial  (convective)  growth  decreases. 

It  was  shown  in  a  previous  QRP  that  either  collisions  or  non-zero 
plasma  electron  temperature  can  stabilize  the  absolute  instabilities. 


v  and  v 
°X  '  0|| 


The  convective  growth  rate  variation  with  03 

.»  parameters  was  studied  In  «R  b.  The  limiting  ca°se  for"  v  app^ach- 
ing  zero  has  been  analyzed  during  this  quarter.  For  small  v"  the 

spatially-growing  backward  wave,  primarily  associated  with  the^'be.m  Is 
given  approximately  by, 


0«  || 


03  -  03  +  i 

c  81/2c 


(7) 


It  is  interesting  to  note  that  when  v  .  0  ,  Eq.  (7)  becomes  identi¬ 
cal  to  Eq.  (b).  This  expression  has  been  checked  with  exact  numerical 
computer  solutions,  and  excellent  agreement  is  obtained  for  m  g  0,8  <n 
lh  Fig.  10.  A,  the  beam  velocl.y  approaches  zero,  the  two  earn  modes  ' 
tend  to  infinity.  The  optimum  biam  velocity  should  be  determined  by 

collisions  or  temperature.  This  quantity  will  be  established  during 
next  quarter. 

It  has  been  po.ntod  out  frequently  that  our  present  theoretical 
e.lort.  are  directed  towards  study  0,  ide.llzed  models  of  be.m/pl.sm. 
interaction,  at id  the  motivation  for  studying  what  becomes  of  the  plane 
whistler  wave  as  plasm,  temperature  and  Inhomogeneity  make  themselves 


-I9- 


(wr/wc) 


625  0650 


•b 


felt  was  briefly  outlined  in  QR  4.1  Although  the  experimental  work 
strives  to  simulate  the  conditions  under  which  the  simple  plane  wave 
theory  is  valid,  the  thermal  spread  and  nonuniformity  of  the  electron 
distribution  will  undoubtedly  affect  the  results.  In  simplest  terms 
it  is  to  be  expected  that  the  desired  whistler  is  itself  modified  there¬ 
to  and  that  additional  wave  types  are  simultaneously  excited  and  coupled 
to  the  whistler  mode.  It  has  accordingly  been  considered  important  to 
develop  the  theoretical  foundations  for  the  analysis  of  wave  propagation 
and  stability  in  a  hot,  magnetized,  inhomogeneous  plasma. 

It  was  explained  in  QR  4  that  this  task  involves  the  combination 
of  the  Boltzmann  equation  for  the  distribution  with  Maxwell's  equations 
for  the  electromagnetic  fields.  A  viable  method  was  expounded  for  ex¬ 
tracting  the  electromagnetic  sources  from  the  perturbations  of  an  equi¬ 
librium  distribution,  by  reformulating  the  system  description  in  terms 
of  inverse  velocity  space.'  The  task  that  remained  was  to  combine  the 
resulting  expression  for  the  sources  in  terms  of  the  fields  with  the 
field  equations  themselves. 

During  this  quarter,  that  task  has  been  carried  out  in  a  formal 
way.  Whereas  it  had  boon  envisaged  that  the  combination  would  result 
in  a  partial  differential  equation  of  appropriately  high  order,  however, 
recognition  has  now  been  given  to  the  fact  that  a  valid  description  of  ' 
the  system  perturbations  in  the  inhomogeneous  case  must  involve  an  in¬ 
tegral  equation,  rather  than  ^he  differential  type.  This  is  a  direct 
result  of  particle  transport  in  the  plasma  and  is  vital  to  the  stability 
problem  because  it  is  the  interaction  of  the  perturbing  wave  with  classes 
of  particles  that  nro  transported  in  resonance  with  the  disturbance  that 
wust  be  expected  to  give  rise  to  wave  growth.  The  resonance  phenomenon 
-hows  up  as  singularities  in  the  integral  elation,  but  it  has  been 
found  possible  to  eliminate  that  clfficulty,  again  by  relying  on  an  in¬ 
verse  velocity  space  description.  The  formalism  developed  for  this  anal¬ 
ysis  has  boon  reported  in  detail  in  a  technical  report  by  P.  Diament 
which  is  at  present  in  the  final  stages  of  preparation.  The  basic  re- 
*  Its  will  be  outlined  here. 
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The  hot,  inhomogeneous  system  is  described  by  its  equilibrium  in¬ 
verse  phase  space  spectrum",  PQ(  6)  ,  which  is  just  the  six-dimensional 
Fourier  transform  of  the  distribution  function;  6  is  the  position  vec¬ 
tor  in  inverse  phase  space.  The  force  fields  to  which  the  particles  are 
subjected,  including  the  applied  magnetic  field  and  the  internal  ambi- 
pola-'  fields,  are  assumed  to  be  describable  by  a  constant  matrix  of  pro 
portionality  to  the  phase,  Y  .  Then  the  Laplace  transform  of  the  time 
development  of  a  weak  disturbance  of  the  spectrum  can  be  shown  to  be 

given  by 


Fi(e,s)  =  il,  J  F0(et  -  60)(et  -  e0)-A6(e0,s  d6eQ/(2*)6 


(8) 


where  L  denotes  Laplace  transformation,  Ag(e,s)  incorporates  the 
perturbing  acceleration  field  and  0t  =  0  exp  Yt  describes  the  unper¬ 
turbed  orbit  in  inverse  phase  space.  This  is  the  solution  of  the  Boltz 
mann  equation  in  terms  of  the  perturbation  Ag  .  The  latter  is  then 
obtainable  by  evaluating  F1(0,s)  and  its  gradient  at  the  origin  in  in 
verse  velocity  space,  A  =  0  .  The  combination  is  expressible  by  the 
pair  of  equations 


a(k,s) 


s2+ 


.2  2 

k  c 


^(k,A,s) 


(9) 


dV 


Fl(*>A>s)  =  J  [F(At'^’s)  +  &>*  X^£o's)]  (at)3  1 


(10) 


where  G^  =  A  X  ^/^A  \  the  arguments  of  FQ  and  G^  are  (k,.  -  At) 

where  6,  =  (k^)  =\a)  exp  Yt  ,  and  «0  and  a  are  essentially  the 
plasma  frequency  and  the  perturbing  electric  field.  This  pair  constitutes 
an  integral  equation  for  the  electric  field,  with  the  plasma  frequency 
as  the  eigenvalue. 

These  results  have  been  specialized  to  the  case  of  a  uniform  mag¬ 
netized  plasma  column  with  a  Maxwellian  distribution  in  velocity  and  a 


-25- 


Gaussian  spatial  profile.  The  final  integral  equation  can  have  either  a 
Gaussian,  or  exponential,  or  Bessel  function  kernel  and  requires  routine 
numerical  analysis  to  extract  the  dispersion  relation  and  stability  cri¬ 
teria.  This  is  planned,  for  the  whistler  mode,  for  the  next  quarter. 
Preliminary  results  already  indicate  that  one  effect  of  inhomogeneity  is 
to  shift  resonances  away  from  the  cyclotron  frequency  by  an  increment 
which  depends  on  the  time  scale  of  traverse  of  a  thermal  particle  across 
a  significant  region  of  inhomogeneity. 

(b)  Experimental  Studies. 

During  the  quarter,  whistler  experiments  have  cot  linued  on  the 
S-band  set-up  using  a  pulsed  reflex  discharge.  Study  of  standing  waves 
indicates  that  the  waves  observed  a  few  hundreds  of  microseconds  after  the 
peak  o*  the  current  pulse  should  be  classified  as  whistlers.  These  waves 
make  smooth  transitions  to  free  space  electromagnetic  waves  at  low  plasma 
densities.  When  plasma  density  is  low,  however,  the  free  space  wave¬ 
length  for  an  L-band  signal  (~  20  cm)  is  too  great  compared  with  the 
plasma  column  cross-section  (7.5  cm  diameter)  to  expect  simple  plane 
wave  theory  to  apply.  Efforts  have  consequently  been  made  to  detect 
whistlers  at  higher  electron  densities.  At  such  densities,  whistlers 
have  been  observed  damped  due  to  collisions.  Quantitative  measurements 
of  their  propagation  have  been  obtained  as  follows: 

The  waves  are  excited  by  short  electric  dipoles  perpendicular  to 
the  magnetic  field  lines.  A  sampling  time  in  the  afterglow  when  the 
plasma  density  is  reasonably  uniform  is  chosen,  and  interf erograms  are 
obtained  by  applying  a  fixed  frequency  and  beating  the  output  signal 
against  a  component  derived  from  this.  Typical  probe  measurements  of 
the  evolution  of  the  plasma  profile,  made  using  a  movable  probe,  are 
shown  in  Fig.  11.  It  is  clear  that  the  density  is  substantially  uniform 
for  times  greater  than  200  ps.  Figure  12  shows  some  interf erograms 
obtained  at  1.2ms.  Since  the  measured  wavelengths  at  ~2.0  GHz  are  of 
the  order  of  1  cm,  plane  wave  theory  should  be  closely  approximated. 

Data  similar  to  that  of  Fig.  12  serve  to  give  the  whistler  disper¬ 
sion  characteristics.  These  are  shown  in  Fig.  13,  and  may  be  compared 
with  plane  wave  theory.  To  do  so,  requires  knowledge  of  the  electron 
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TIME  IN  THE  AFTERGLOW! 


■2', 


FIG.  11  Pulsed  reflex  discharge:  Variation  of  electron  density 


density.  This  was  obtained  by  calibrating  the  Langmuir  probe  measure¬ 
ments  by  means  of  a  K-band  interferometer.  Good  ngroomont  is  indicated 
in  the  figure  for  the  range  covered  by  the  mea suromonts .  To  oxtend  them 
to  lower  (0)/^)  ,  highor  magnetic  fields  are  needed.  Furthor  measure¬ 
ments  will  be  made  next  quarter  in  the  new  magnetic  field  system, which 
gives  cyclotron  frequencies  up  to  K-band.  The  status  of  this  set-up  is 
as  follows.  The  magnet  system  that  has  delayed  this  project  has  boen 
tested  and  is  now  available.  It  is  shown  in  Fig.  14.  It  con  provide  a 
steady  field  of  up  to  7.7  kGauss  which  is  uniform  to  ±  0.25^  over  a  dis¬ 
tance  of  40  cm.  The  maximum  AG  ripple  of  the  magnetic  field  is  less 

than  1$,  and  the  ripple  at  full  power  is  less  than  0.1<£.  The  magnet  dis 

sipates  nearly  150  kW  at  maximum  power. 

A  vacuum  system  (not  shown  in  figure)  has  been  assembled,  and  exper¬ 
iments  now  under  vay  in  the  set-up  are  directed  towards  obtaining  a  high 
density  rf  discharge  at  a  few  MHz.  The  rf  transmitter  that  will  be 
used  to  produce  the  plasma  is  in  the  final  testing  stages.  The  various 
probes  and  interferometers  which  will  be  used  as  diagnostic  devices  hove 

already  been  constructed.  The  discharge  will  be  produced  by  the  rf 

transmitter,  operated  in  a  pulsed  mode,  and  allowed  to  decay.  A  maxi- 

13  3 

mum  plasma  density  of  approximately  10  /cm  is  expected. 
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FIG.  14  K-band  magnetic  field  system  for  use  in  whistler 
amplification  studies. 
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V 


FUTURE  PROGRAM 


Most  of  the  details  of  our  program  for  the  coming  quarter  have  been 
dealt  with  in  the  relevant  theoretical  and  experimental  sub-sections  of 
Sections  II-IV.  Summarizing,  the  program  is  as  follows: 

(i)  Beam/plasma  amplification  with  transverse  modulation  — 
Theoretical  work  will  continue  on  beam/surface  wave 
interactions, first  for  the  m  =  0  mode,  to  check  the 
experimental  results  already  obtained,  then  for  the 

m  =  1  mode.  It  is  anticipated  that  our  measurements 
on  the  sealed-off  tubes  constructed  so  far  will  be 
completed,  and  that  further  studies  will  be  made  in  a 
more  flexible  continuously-pumped  system. 

(ii)  Electrostatic  wave  amplification  in  magnetoplasmas  — 

Further  measurements  on  the  noise  spectrum  due  to 
magnetoplasma  wave  excitation  by  electrons  with  trans¬ 
verse  energy  will  be  made,  first  with  a  view  to  iden¬ 
tifying  the  various  frequencies  so  far  observed,  then 
with  the  aim  of  verifying  the  theory  quantitatively 

for  a  delta-function  beam  interacting  with  a  cold  plasma. 

(iii)  Electromagnetic  wave  amplification  in  magnetoplasmas  — 

Our  studies  of  the  relevant  dispersion  relations  will  be 
extended  numerically.  It  is  hoped  that  the  propagation 
measurements  will  be  completed  in  the  S-band  magnetic 
field  system.  They  will  then  be  extended  in  the  K-band 
system  preparatory  to  introducing  an  electron  beam  with  trans¬ 
verse  energy  to  excite  wave  growth. 
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